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Abstract 

Radiative heat exchange at the nanoscale presents a challenge for several areas due to its scope and 
nature. Here, we provide a thermokinetic description of microscalc radiative energy transfer including 
phonon-photon coupling manifested through a non-Dcbye relaxation behavior. We show that a lognormal- 
like distribution of modes of relaxation accounts for this non-Debye relaxation behavior leading to the 
thermal conductance. We also discuss the validity of the fluctuation-dissipation theorem. The general 
expression for the thermal conductance we obtain fits existing experimental results with remarkable 
accuracy. Accordingly, our approach offers an overall explanation of radiative energy transfer through 
micromctric gaps regardless of geometrical configurations and distances. 

Introduction 

There is a general consensus that the search for clean sources of energy with no climatic and environmental 
impact constitutes a major strategic objective at present. In this sense, nanoscale thermal radiation 
conversion offers a source for intensive clean energy generation. Thermal radiation has always been an 
active field of study. Throughout the 19 th century, great scientists (Boltzmann, Stefan, Rayleigh, etc.) 
dedicated a considerable effort to this problem which was completely solved at the turn of the 20 th 
century due to Planck's contribution to the founding of quantum mechanics. After this, it seemed to be a 
well-established fact that the maximum power extracted from a hot body depended on the temperature 
as T 4 . 

However, recently, thanks to modern technological advances, it has been shown that energy exchange 
through thermal radiation at nanometric distances breaks by several orders of magnitude the limits 
posed by Stefan-Boltzmann law for black body radiation. Moreover, near-field thermal radiation is 
approximately monochromatic and reveals itself coherent in space and time, which may lead to stationary 
interference phenomena in a micro-cavity. Therefore, this monochromaticity and coherence along with 
the overcoming of Stefan-Boltzmann limits, all of these distinguishing features together confer near-field 
radiation a great potential for future applications in nanotechnology and, as we have said at the beginning, 
energy conversion as well. Several reviews on this issue have been written recently, see Ref. [IH3] by way 
of example. 

An ever-increasing number of investigators has marked the recent history of the research on near- 
field radiation. Polder and van Hove [3] first studied heat transfer between two objects at nanometric 
scales maintained at different temperatures by following a stochastic or fluctuational electrodynamics 
formalism established by Rytov et al. [5]. Recently, it has been emphasized [6] that surfaces modes as 
included in the solution of Maxwell's equations in Ref [4] can greatly enhance the heat flow. Experiments 
showing that the heat flow at the nanoscale is indeed greater than the blackbody radiation limit among 
materials supporting surface modes were reported in Ref [7J[5] and between two gold surfaces in Ref [5]. 



2 



Likewise, Pcndry [lOj gave a simple derivation of the expression found by Polder and van Hove and 
some interpretation in terms of heat transfer channels in addition to a discussion of the maximum heat 
flux. The problem was reformulated by using a Landauer-likc approach by Biehs et al. Finally, 
Sasihithlu and Narayanaswamy |12| performed a discussion of the proximity approximation. All these 
approaches have something in common: the linear response regime and, consequently, the fluctuation- 
dissipation theorem (FDT), whose validity is not guaranteed at this level. Thus, an approach based on 
the contributions of fluctuating dipole effects seems to be the heart of a considerably simplified treatment 
of energy transfer at the nanoscale [7]l8]ll3j. These approaches include a wide range of phenomena in 
which the energy between molecules is dominated by dipole-dipole interactions, also known as the Forstcr 
energy transfer [14]. 

Nonetheless, as two nanostructures thcrmalizcd at different temperatures come closer to each other, 
the distribution of charges and currents becomes asymmetric and therefore, defies description in terms 
of dipolar interactions. Hence, it becomes clear that one must bear in mind higher order effects beyond 
the dipole |15j and also include other contributions to thermal conductance quite common in disordered 
amorphous materials, leading to a generalization of the FDT mentioned above [16] . 

In this article, we will go deeper into these aspects by showing how phonon-photon coupling effects 
account for microscale radiative energy transfer by considering non-Debye relaxation due to the excess of 
modes in the low and high frequencies in the bulk material |15H17j . In the current literature on disordered 
systems, the non-Debye refers to an excess of modes of vibration over the Debye level observed in inelastic 
light (Raman) and neutron scattering. From here, we obtain a general expression for the heat transfer 
coefficient including both Debye and non-Debye contributions, providing an overall explanation of the 
energy transfer through micrometric gaps. The findings of our theory fit existing experimental results 
with a high degree of accuracy. 



Our theoretical framework has been described in previous publications |16pi7] and is briefly summarized 
here. We consider a gas of quanta distributed in phase space according to the probability density p(T, t), 
where T = (p,x), p = |p|, and x = |x|, with p and x being the momentum and position of a quanta, 
respectively. Here, we must pay attention to the fact that p(T,t) possesses dimensions of h~ 3 , h being 
Planck's constant. The thermodynamic description that we propose entails the formulation of the second 
law of thermodynamics, which can be carried out by means of the Gibbs entropy postulate [18] . 



This equation gives us the nonequilibrium entropy of the gas of quanta plus the bath, with S eq being the 
equilibrium entropy and p eq the equilibrium probability density. 

In general, entropy is produced due to irreversible processes, in such a manner that irreversible 
processes in nonequilibrium systems are described by means of currents, thermodynamic forces (affinities), 
and the entropy production rate, which is always positive. It is precisely this positive character of 
the entropy production that enables us to derive relaxation equations. Since the probability density is 
conserved the existence of a density gradient dp(T,t)/dT = (dp/ dx, dp/ dp) yields a current J(T,t) = 
(Jx,Jp) which unleashes a relaxation process dp/dt = —dJ/dF. In addition, this current satisfies the 
relation J (T,t) = — D(p) • dp(T,t)/dT, derived from the entropy production, which can be obtained 
from Eq. ([I]) 16, 18 . Here, D(p) is a material-dependent quantity, the matrix of diffusion coefficients, 
satisfying Onsager's symmetry principle. In general, due to the tensorial character of the diffusion matrix, 
both currents J x and J p are coupled 



Methods 




(i) 



J a = -D a p — ; with a,/3 = x,p, 



(2) 
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where D a @ are the diffusion matrix components. 

The main contribution of this article is to apply this formalism to the description of the radiative heat 
transfer between a nanosphere and a plate at different temperatures, T\ (hot) and Ti (cold), separated 
by a distance d (see Fig. [TJ. We assume that heat transfer results from two different mechanisms. (I) On 
one hand, we consider a conventional radiative heat exchange involving the dynamics of quasiparticles as 
the result of two simultaneous processes: elastic emission and absorption of hot photons from the medium 
at T\ and elastic emission and absorption of cold photons from the medium at these processes also 
involve the presence of surface modes [6j[T0]. (II) The second mechanism behind heat transfer has 
to do with the excitation of coupled resonant modes from the collection of acoustic states related to 
defective soft structures in a disordered regime [T9H2T] . This process encompasses inelastic scattering of 
the impinging radiation, which is linked to nonequilibrium contributions. In both scenarios, there is no 
diffusion in configuration space since quanta are massless particles [16j . and as a consequence J x = 0, 
which brings about J(T,t) — > J p (T,t). Hence, from Eq. ^ one obtains the appropriate diffusion 
current in momentum space 

D PP D XX \ dp{T,t) _ h d P (T,t) 



where r(p) is the relaxation time, which depends on the diffusion matrix components. 
Near-field analysis 

In the near-field regime, confinement of the electromagnetic waves in a micromctric gap separating neigh- 
boring nanostructures introduces peculiar effects in the spectrum of the thermal radiation. Here is where 
the collective modes (phonons) excited in the material by the impinging radiation come into play. To 
discern whether confinement effects are important or absent we shall assume a cut-off wavelength, the 
thermal wavelength of a photon At = hc/ksT, which is proportional to the Wien's displacement law 
through a proportionality constant, i.e. At = P/T, where here (3 as 2.82143937212. Actually, when 
d 3> At we have a blackbody spectrum of radiation. One may wonder what happens when d < At- 
Since according to Heisenberg's principle Ax Ap ^ h, assuming that the maximum value of Aa; is d, one 
obtains Ap 2i h/d, and thus the minimum value of Ap is h/d, leading to AE ^ hec/d, with < e < 1. 
Here, e is, precisely, the inverse of the refractive index. Hence, not all frequencies are possible and the 
frequency of resonance ujr — 2nec/d appears. 

In these circumstances, integrating by parts the resultant continuity equation gives us 

|)*-* M 

with J p r~j p/th 3 , where we have assumed that J p (T, t) — at p = — oo. Performing a second integration 
of Eq. ([3]) through from 1 to 2, we find the net current 

J (P,t) = T77t(P2-Pi), (5) 

T*(t) 

where pj(p, t) = p (j>, x = Xj, t) is related to the population of quasiparticles at Xj and pi + p2 = 1- More- 
over, J(p,t) = (1/t*) Jj 2 r(p)J p dx, with r*(i) = j pr(p)dpdx corresponding to a hierarchy of relaxation 
times ubiquitous in complex systems. 

In the stationary state, once the system has thermalizcd at temperatures T\ and T2, pj(p,t) — > 
2N(uj,Tj)/h 3 ; the factor 2 comes from the polarization of a photon and N(uj,T) being the averaged 
number of quasiparticles in a elementary cell of volume ft 3 of the phase-space given by Planck's distri- 
bution [52], N(w,T) = 1/ [exp(Sw/fcT) - 1]. Besides, r*(t) — > t*(lu); therefore leading to a stationary 
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value, i.e. J(p, t) — > J st (w). Thus, the heat flow Q can be obtained from the sum of all the contributions 



as 

Q = ec.J st (io)dp = -h 3 ec / J st (uj)g (w) du, (6) 

where p = (huj/sc)^lp, with Jl p being the unit vector in the direction of p, and the distribution of 
frequencies is given by 

9{UJ) = ^cf- (7) 

At this point, it is worthwhile to make a short digression about the physical meaning of the time 
scale t* (lj). It has been known for a long time [23] that for most condensed systems in time-dependent 
fields, the orientation polarization behavior can, as a good approximation, be characterized by a relaxation 
time distribution (r* (i)); this behavior is generally meant as dielectric relaxation. In harmonic fields, this 
implies that the complex dielectric permittivity in the frequency range corresponding to the characteristic 
times for the molecular reorientation can be written as 

poo 

e(u) = e 00 + (e s -e 00 ) ^ / dr, (8) 

J 1 + IUJT 

where e s is the static dielectric constant and is the permittivity at the infinite frequency. Here, 
the relaxation time distribution Q it) satisfies the normalization condition, J °° £ (r) dr ~ 1 and Eq. ([S]) 
constitutes a generalization of the Debyc treatment based on Clausius-Mossotti equation [2"4] 

e M = £oo + (e s - £oo) — r~. — • (9) 

1 + ILOT 

This Debye's equation, Eq. ©, follows after Fourier-transforming the relaxation function (j> (i) = 
exp (— t/r) which coincides with the normalized dielectric function (e s — e(t)) / (e s — e^). Note that 
if we assume a time-independent time scale r in Eq. ([5]), after integration we shall obtain the Debye 
relaxation function in terms of exp (— t/r). Hence, we can understand t* (lu) as defined through the 
relation 

rj^L dTS — l_ (10) 

which shows that the Callen-Welton FDT [55] is not valid at this level. In fact, unlike here, the FDT is 
related to decaying equilibrium fluctuations characterized, precisely, by a single relaxation time. 

Now, let's return to the main topic after that brief digression. Note that according to Eq. (|6]), in the 
limit d — > co, 

1 f°° 

Q=-h 3 ec J st (uj)g(uj)du, (11) 
1 Jo 

giving us the blackbody radiation limit provided t*(oj) oc lo, Q cx (fc B Ti) 4 - (k B T 2 ) ■ On the other 

hand, in the limit d — > 0, Q — > which in contrast to the descriptions based on evanescent surface waves 
avoids divergences in the heat flux in a self-consistent way. Nonetheless, for finite d we can rephrase the 
expression of the heat current given by Eq. ([6]) introducing a new varible uj = 1/s 



1 /""« ds 

Q = ^h 3 ec J st [u(s)]g[uj(s)}^ 

L Jo s 

h 3 ecujR 



Jst [w (x^r 1 )] 9 [w (x«>r)] , (12) 
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where the mean value theorem has been used to approximate the integral, with < x < 1. Hence, since 
w (xw^, 1 ) = x^ 1 ^R, Eq. (fT2|) reduces to 



Q = 



h 3 ecwn 



Jst (x 1 ur) 9 (x 1 ur) ■ 



In terms of Planck's distribution, Eq. (|13[) can be rewritten 

hecoj R g (x~ 1 ur 



Q = 



X 2 t* (x 



[N (x-^Tx) - Nix-^Ti)] , 



(13) 



(14) 



while t*(lu) must be determined in a more general way, by scrutinizing the interaction processes between 
light and bulk material. 

For first order in the temperature difference AT, 



k B ATec g {x ur) 
X t* (x^^r) 



hX 1 ujR/2k B T Q 



smh(Hx 1 ujR/2k B T ) 



(15) 



with AT = Ti — T2, and Tq = (Ti + T2) /2. Thus, it follows that the heat transfer coefficient, the quantity 
usually measured in experiments and defined by Qj AT, is given through 



H(d,T ) 



X 1 ur 



Tr 2 e 2 c 2 x t* (x 1 Ur) 



hx-^R^kBTp 
sinh (hx~ 1 WR/2k B T Q ) 



(16) 



When the mechanism of heat exchange is through elastic collisions, which is similar to Rayleigh 
scattering, it is known that the intensity of radiation is proportional to lu 4 |26j . Therefore, 



T*{^)- l =T- l /A, 



(17) 



where the time scale r is a material-dependent parameter. On the other hand, regarding the inelastic 
contribution to the near-field heat exchange, this is the analogue to the Raman scattering of light. In 
this case, the distribution of modes presents anomalies which result from states located at a lower energy 
region [31]. The Raman spectra is fitted using a lognormal function first proposed by Denisov and 
Rylev [19]. This lognormal distribution is a statistical model, which can describe collective motions 
causing extremely slow structural relaxation, thereby fitting the non-Debye anomalies |27j . Therefore, 
this accounts for the high nonlinear behavior of the thermal conductance between both materials. For the 
proposed case, we assume that the density of vibrational states is achieved through the use of a lognormal 
distribution, which corresponds to 



4v27rcra;o 



exp 



In 2 (u/u ) 
2a 2 



(18) 



Here, ujq (characteristic frequency) and a (standard deviation) are two fitting parameters characterizing 
the lognormal distribution. The lognormal in Eq. (|18p stems from the existence of a hierarchy of relaxation 
mechanisms in the material, related to the presence of collective effects. This distribution, results from 
the fact that the energy of the system consists of a large number of contributions and the application 
of the central limit theorem of probability theory (see Appendix SI,). In view of the properties of the 
lognormal distribution, it must be noticed that 



1 



-00 WT*(W) 



d(ln (w/wq)) 



1 



4t o w ' 



(19) 



being a kind of closure relation for the relaxation times. 



G 



Hence, the excitation of a mode constitutes a photoinduccd cooperative phenomenon. It is plausible 
to assume that the cumulative effect of incident photons ends perturbing the material, thus triggering 
collective oscillations. In addition, the lognormal accounts for an excess of ways of adsorbing energy by 
the system with respect to the ways obtained when merely the Debye squared-frcquency law describes 
the relaxation. Therefore, this provides a reasonable description of the non-Debye law [20][2Tj, which as 
in Raman scattering also becomes manifests in radiation problems. 

Accordingly, the heat transfer coefficient, Eq. (|16j) . results from the addition of the elastic and inelastic 
contributions mentioned above 

if (d, To) = 4 

T o 

where t' = % 3 r and v = ex 1 - In Eq. ([20]) , the first term inside the square brackets corresponds to the 
usual contributions found up to now in the current literature [6], taking into account surface phonon- 
polaritons owing to the presence of evanescent waves close to the interface. As we have mentioned above, 
the second term takes into account cooperative phenomena, becoming manifest through the existence 
of collective modes of vibration in the system, which appear in the density of states. Consequently, a 
more general formulation for the radiative heat transfer problems must come from the superposition of 
Debye and non-Debye relaxation mechanisms, combining in this way the contributions from the material 
surfaces as well as the bulk. 



(27T) 



1/2, 



crujQd 



■ exp ■ 



In (2-kvc/ 'loqcI) 



( hvc/2kBT{jd 
Vsinh {hvc/2k B T {) d) 



(20) 



Results and discussion 

The thermal conductance is obtained by integrating the heat transfer coefficient over the surface of a 
sphere of radius R divided by its area {A = 47ri? 2 ), i.e. G = (l/A) J" H[d{r),T ]2Trrdr. Since the 
distance d between both surfaces depends on geometrical characteristics, the local distance between the 
sphere and the plane surface must be measured through the local radius r, assuming the effective distance 
as d(r) = d + b + R — V R 2 — r 2 , with b being a surface roughness parameter [7] . 

We have calculated numerically the surface average by adjusting the parameters mentioned above to 
the experimental results obtained in the Ref [5] , which takes into account only the near-field contribution 
when decreasing the sphere-plate distance. In Fig. [2] we show the near-field conductance fitting the 
values of the integral of Eq. (|20[) to the data for glass-glass [8] and gold-gold materials [9]. For both 
the material we have used (r^)" 1 = 2.1 x 10 7 Hz, w = 1.7 x 10~ 13 Hz, and a = 6.0 x 10 18 . For the 
glass-glass material we have obtained v = 3.2 x 10~ 3 and for gold-gold v = 7.9 x 10~ 4 . The near-field 
heat transfer described using evanescent waves as solutions of classical electrodynamics equations leads 
to heat flux divergences as the gap vanishes |28| . However, in our approach this divergence does not 
occur, reaching a constant-conductance value as the distance between the nanoparticles decreases (shown 
in the inset). In essence, this is due to the fact that the time relaxation distribution herein described 
by lognormal distribution incorporates two effects: (I) Phonon branches in a real structure affecting the 
density of states in different frequency regions, similar to actual behavior observed in metals [29) : and (II) 
the vibrational modes in the bulk material absorbing the energy excess. Hence, since our assumption of a 
distribution of relaxation times accurately describes the dynamics of radiative systems at the microscalc, 
we conclude that the FDT in the Callcn-Wclton formulation [25] is not applicable at the nanoscale and 
must be modified. By going beyond the Debye theory, a way for this generalization is offered here. 

Conclusions 

In summary, we have evaluated thermal conductance in the near-field, giving a thermokinetic description 
of some experiments involving heat radiation through a very narrow gap. Although near-field radiative 
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transfer is a highly complex phenomenon, we have been able to provide a unified and highly accurate 
explanation of heat exchange processes at the nanoscale. Our theory covers all distances from the far-field 
up to contact. Since the experiments examined may involve a great variety of nanostructures, our theory 
possesses a wide scope of applications. The general methodology presented here may also be used in the 
study of other heat exchange processes such as those occurring in phonon systems and in the analysis of 
thermal contributions to Casimir forces, even in charge conduction problems in nanosystems. 
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Figure Legends 



Figure 1. Schematic diagram of the radiation exchanged via conventional radiative 
transfer. Elastic collision of photons with atoms or molecules of materials and phonon-photon coupling 
contributions between a sphere and a plate maintained at different temperatures, T\ and T2, separated 
by a distance d. 



225501. 
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Figure 2. Sphere-plate near-field heat transfer coefficients between a gold (or glass) 
sphere and a gold (or glass) substrate versus gap distances. The data are from Ref [9] for the 

50 diameter spheres. The dotted lines are comparisons with the theoretical predictions from the 
proximity theorem. The inset shows a non-divergent regime as the gap vanishes. 



Appendix SI: Lognormal distribution 

This Appendix is devoted to the derivation of the expression of the lognormal distribution corresponding 
to Eq. (18). The energy of the system results from the combined effect of a large number of independent 
inputs e Q , ei + ... + e„. Let us assume that the effect on the energy of the input e v is proportional to e v 
and to the cumulative effect E v of the v — 1 previous inputs, E v +\ = E v + e v E v . Whence 



e 1 + ... +e n = 



E, 



E„ 



E„ 



de 



1 ^ E 
— log — 

k 5 E l 



(21) 



For large n, the distribution of the sum is given by a Gaussian, according to the central limit theorem 

2 



p (u>) ~ cxp 



1 w 
-log — 

k U!l 



/2a 2 



(22) 



where E = hto. Equivalcntly 

p (uj) duj 



ku>i 



2lT(TUJ 



exp 



1 u 
- log — 

k u)\ 



/2a" 



du. 



(23) 



This distribution depends on two empirical parameters a and u>±. 



